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^ ! Abstract 

■ The production of various heavy-quark-containing baryons in electron-positron annihilation is 

considered. On the basis of exact formulas that we obtained previously within full perturbation 
theory, new numerical calculations of the respective cross sections are performed, and simple ap- 
J> I proximate expressions are then constructed for the results of these calculations. The dependence of 

the total cross sections on the masses of constituent quarks is discussed. The application of the Pe- 
terson fragmentation function and a Reggeon-type fragmentation function to describing differential 
^ . cross sections is analyzed. 

O 

1. Introduction 

Qh! Investigation of the mechanisms responsible for the production of hadrons containing 

I heavy quarks is of interest from the theoretical point of view since this provides the possibility 
for further testing QCD - more precisely, our understanding of it. In this way, one tests both 
its perturbative aspect used to describe the simultaneous production of several quark pairs 
and nonperturbative models constructed on the basis of QCD for bound states. We recall 
that, even in cases that are the simplest at first glance, the results of calculations appear 
to be in an unexpected contradiction with experimental data, as was, for example, in the 
hadronic production of J/ip particles. At the same time, derivation of theoretical estimates 
for relevant cross sections is of importance for practical purposes, such as those associated 
with planning searches for such particles and investigations of their properties. 

Available calculations of the cross sections for the production of baryons containing heavy 
c and b quarks rely, as a rule, on considering the production of respective diquarks, this 
corresponding to the fourth order of perturbation theory. A detailed review of the results 
obtained in this way and an exhaustive list of relevant references can be found in [1] . In the 
sixth order of perturbation theory [(9(a;^ag)] , the total and differential cross sections for 
the production of multiply heavy baryons flscb and flccc at the Z pole in electron-positron 
collisions were calculated in our previous studies [2, 3]. For the squares of relevant matrix 
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elements, wc obtained exact analytic expressions, which, as might have been expected, are 
very cumbersome (they exist only in the form of computer codes). As a result, numerical 
calculations with these matrix elements would be extremely time-consuming. 

By using these expressions and performing a series of new numerical calculations of vari- 
ous cross sections for the production of baryons containing three nonidentical quarks, we try 
here to impart, to the emerging results, a broader content, simplicity, and adaptability in the 
possible future application to constructing estimates for planning experiments. Specifically, 
we find, first of all, a simple approximate dependence of the total cross section for baryon 
production in electron-positron collisions on the mass of each constituent quark. On the 
basis of the concept of fragmentation, we then approximate the differential cross sections 
with the aid of the Peterson function for various sets of quark masses. 

In this study, we also analyze some aspects of the description of the differential cross 
sections for the production of flscb and Qccc baryons in terms of a Reggeon-type fragmentation 
function. 



2. Dependence of the total cross sections for baryon production on 
consistent- quark masses 

Let us consider some properties of the production of a qiq2q3 baryon consisting of three 
nonidentical quarks qi, q2 and ^3 and having a mass M, a momentum p, and an energy E 
at the Z pole in electron-positron collisions. We assume that the quark masses mi, 1712 and 
rris differ markedly from each other; that is, 

ml <^ ml <C ml- (1) 

In QCD, the elementary process corresponding to the production of such a baryon is 

e+ + e" ^ gi(pi) + 52(^2) + q^iPs) + Qiipd + ^2(^5) + QsiPe), (2) 

where the quark and antiquark 4-momenta are indicated in parentheses. 

The formation of the baryon from three quarks is described in the well-known nonrel- 
ativistic approximation [4-6], whose details required for our purposes are as follows. First, 
the velocities of the quarks forming the baryon are assumed to be identical. Second, the dif- 
ferential cross section for baryon production is obtained from the standard differential cross 
section for the process in (2) by replacing the phase space of three quarks by an expression 
proportional to the baryon phase space; that is, 

d'pi d^P2 d'ps ^ \m\' 

{2tt)^ ■ 2Ei {2n)^ ■ 2E2 {2nf ■ 2E3 ^ {2tt)^ ■ 2E' 

where ip{0) is the value of the wave function at the point where relative coordinates of all 
three quarks are zero. 

In each of the Feynman diagrams corresponding to the process in (2), one can indicate 
a virtual gluon g such that it transforms into a quark-antiquark pair qi{pi)qi{pi+3) without 
emitting a new gluon g'. The denominator of the propagator of the gluon g - it has the form 
{Pi +Pi+3y ^ attains a minimum of 4m^ at = Pi . But if a virtual gluon g transforms into 



2 



two quark- ant iquark pairs qi{Pi)qj{pj)(ji{Pi+3)(jj{pj+3), the denominator of its propagator has 
a minimum of 4(mj + rrij)^. Taking into account the inequalities in (1), we can deduce from 
the above that the leading contribution to the amphtude of the process in (2) comes from 
the diagrams where the production of quark- antiquark pairs proceeds hierarchically from 
the heaviest to the lightest. This sequence that ends up in the formation of a baryon from 
the quarks qi , q2 and q^ is generally referred to as the fragmentation of the quark q^ into a 
?i?2?3 baryon. It is often described analytically in the form 

ag,q,- D{z), (4) 

where (Jq^q^ is the total cross section for the process e+e" — > q^q^ and D{z) is the respective 
fragmentation function. For the variable z, one usually takes the quantity Xp — p/pmax or 

Xe = E / -Emax- 

It is obvious that the amplitude of any process is a homogeneous function of the 4- 
momenta P^ , the masses Mj , and the widths Vj of real and virtual particles involved in the 
process. Therefore, the total cross section or one differential cross section or another can be 
represented in the form of the product of some power of the total energy ^/s and a function 
of the reduced 4-momenta Pj/ y^, the reduced masses Mj/ a/s, and the reduced widths Tj/ s. 
From here, it follows, among other things, that the fragmentation function D[z) appearing 
in (4) depends parametrically on the reduced masses mij ^/s of the product quarks. In the 
following, we will write the reduced masses explicitly only in the logarithmic factors on the 
right-hand side of formula (6) (see below). 

In comparing experimental results obtained in electron-positron collisions with some 
fragmentation function, attention is given primarily to its form depending on one or two 
parameters but not to its normalization. 

For want of experimentally observed events involving the production of qiq^qa baryons 
in electron-positron annihilation, it is reasonable to focus on the total cross section-namely, 
on the dependence of the total cross section on the masses of the quarks qi, q2 and q^ . The 
choice of a simple algebraic expression representing this dependence is based on the fact 
that the square of the matrix element of the process being considered is similar, in some 
respects, to a rational function of the momenta of product particles, its denominator at the 
minimum involving constituent quark masses and their sums as factors. It is well known 
that the integral of such a function can generally include logarithmic terms. 

By using relation (3) and the results of numerical calculations of the total cross section 
for the production of q-^q^qz baryons at the Z pole in electron-positron collisions for six sets 
of masses mi,m2 and rris for the same set of electroweak-interaction coupling constants (such 
as that for the production of scb baryons), we arrive at the formulas 

n l^(0)P r 

(mi 777.2 + "73)"= 



C 

m\m2 \4mi J \fn3 
where ^=91.2 GeV and C = (0.0407 ± 0.0006) pb. 
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In Table 1, we present the sets of masses mi, m2 and m^; the values of G that are 
obtained from a Monte Carlo calculation of the integral of the square of the matrix element 
for the process in (2) over the phase space of four final particles (QCD column in the table); 
and the values of G that are obtained by formula (6). 

We believe that the use of the factor Zn(\/i/(4mi)) in formula (6) is quite justified 
empirically. It seems plausible that the quantity G depends only slightly on the mass m,^ 
of the heaviest quark; relying on the results of our calculations exclusively, we cannot be 
confident, however, that the G depends on rris through the factor ln{^/m'i), as follows from 
Eq. (6). 

Table 1 



mi 

(GeV) 


m2 

(GeV) 


mz 

(GeV) 


QCD, G 
(pb ■ GeV^) 


Formula (6), 
(pb ■ GeV^) 


G 


Parametr e 
from (8) 


0.5 


1.5 


4.8 


(7.85 ±0.16) 


10-1 


(8.15 ±0.14) 


•10-1 


0.124 ±0.016 


0.3 


1.5 


4.8 


(2.58 ±0.07) 


10° 


(2.57 ±0.04) 


•10° 


0.098 ±0.012 


0.075 


1.5 


4.8 


(5.59 ±0.20) 


10^ 


(5.42 ± 0.09) 


•101 


0.066 ± 0.008 


0.01 


1.5 


4.8 


(4.29 ±0.12) 


10=^ 


(4.12 ±0.07) 


•10=^ 


0.048 ± 0.006 


0.01 


0.5 


4.8 


(3.90 ±0.28) 


10^ 


(3.71 ±0.06) 


• 10^ 


0.016 ±0.002 


0.01 


0.5 


1.5 


(5.10 ±0.53) 


10^ 


(5.18 ±0.09) 


• 10^ 


0.24 ±0.08 



Strictly speaking, the dependence of the total cross section on the quark masses mi, 
m2 and m^ is not exhausted by the explicit expressions in formulas (5) and (6), since the 
baryon-state wave function at the origin, '0(0), must change in response to a change in the 
masses. The dependence of ^(0) on the masses of constituent quarks is determined within 
potential quark models, which are not discussed in the present study. On the basis of the 
numerical values of |'?/'(0)p that are presented in [7] for six sets of quarks gi, g2, and ^3, we 
can nevertheless estimate cross sections by using the approximate expression 

|V'(0)p -Dm2^m3, if i/ mi < m2 < ma, (7) 

where D = 0.065 • lO'^ GeV^ ^ 

Taken together, relations (5)-(7) indicate that the total cross section for the production of 
Q1Q2Q3 baryons in electron-positron annihilation is highly sensitive to the mass of the lightest 
of three quarks. This circumstance is especially important in the case where, for the quark 
qi, one takes a -u or a quark, since, from the point of view of simple nonrelativistic concepts, 
their masses can be varied within rather wide intervals, from 50 MeV (in pions) to 300 McV 
(in nucleons). In order to estimate cross sections for the production of baryons containing 
two heavy quarks and au or a d quark, we set m^ = m^ = 100 MeV in Eqs. (5)- (7), bearing 
in mind that, according to the approximation specified by Eq. (6), these cross section are 
determined to within a factor of 10. 

For want of something better and, to some extent, as a continuation of the strategies 
adopted in [7], we propose extending the procedure of the present study and of previous 
studies [2, 3] to the case of deriving estimates for the production of baryons containing 
one heavy quark c or 6 , a strange quark s, and a light quark u or d-namely, we propose 
supplementing the perturbative part of calculations (sixth order of perturbation theory) with 
the following nonrelativistic nonperturbative elements: the assumption of equal velocities of 
the quarks fusing into the baryon in question, relation (3), and the approximation specified by 
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Eq. (7). What we inherit from nonrelativistic potential models reduces to the extrapolations 
in (3) and (7). It is reasonable to indicate here that, in contrast to u and d quarks, a 
strange quark of mass = 500 MeV leads to an acceptable level of agreement with naive 
nonrelativistic expectations for the masses of the meson and baryon ground states [ss (0) 
and sss (Q). respectively]. 

All estimates that can be obtained on the basis of formulas (5)-(7) by setting rrib = 4.8 
GeV, rUs = 0.5 GeV, and m„ = = 0.1 GeV are given in Table 2 (the symbol q in the 
subscripts there stands for a or a c? quark). The factor C in (6) is set to 0.0407 pb in 
calculating cross sections for the production of ^Igcb , '^qcb , and Eqsb baryons and to C = 
0.0407 pb -{{g^f + + 9aY = 0.0317 pb in dealing with Eg^c baryons. 

Table 2 



Baryon 


103- 1^^(0)12, GeV^ 


o'tot, fb 


^scb 


0.57 


0.0097 


^qcb 


0.57 


0.40 


^qsb 


0.11 


0.98 


^qsc 


0.034 


2.2 



We would like to bring to the attention of the reader that, in [8], the quantity |^(0)p 
was calculated at the following values of the constituent quark masses: = 5.29 GeV, rric 
— 1.905 GeV, rUs = 0.6 GeV, and m„ = rrid = 0.3 GeV. In our previous studies [2, 3], we 
employed the values of |-?/'(0)p from [8] without introducing any corrections. In the present 
study, all cross-section values, including those that are given in the figures, are calculated 
by using Eq. (7). 

It should be emphasized that the approximate formula (6) was obtained for the con- 
stituent quark masses obeying the inequalities in (1). But if these inequalities are not 
satisfied and if, in addition, there are quarks identical in fiavor among the quarks in ques- 
tion, formula (6) cannot be used even for a rough estimate of cross sections. Indeed, it 
can be seen that, for the production of an flccc baryon in electron-positron annihilation, a 
direct numerical calculation of the quantity G on the basis of formula (5) gives the value of 
2.27 pb GeV^ while expression (6) yields 0.090 pb GeV^ The reasons for so significant a 
discrepancy between the values of G are rather obvious: if the inequalities in (1) hold, the 
main contribution to the cross section comes only from a few Feynman diagrams, but, if all 
three masses mi, 1712, and ms are close to one another, all diagrams make approximately 
equal contributions (for the production of an Qccc baryon, there are 504 such diagrams); 
additionally, interference effects arise if there are identical quarks. 



3. Differential cross section and Peterson fragmentation function 

Let us now proceed to consider a simple algebraic description of the differential cross 
sections for the production of gig2?3 baryons in electron-positron collisions on the basis of 
the fragmentation approach. Experimental data on the production of heavy-quark hadrons 
in electronpositron annihilation are usually approximated in terms of the Peterson function 

[91 ^ 
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where e is a phenomenological parameter. We used this function in [2] to describe approxi- 
mately the production of flscb baryons in electron-positron collisions. 

We performed complete numerical calculations of the differential cross sections for the 
production of qiq2qz baryons for six sets of constituent quark masses (see Table 1) and then 
determined the values of the parameter e of the Peterson function (8) that ensure the best 
agreement between the form of this fragmentation function and the form of the differential 
cross sections da/dxp. The resulting values of the parameter are given in the last column of 
Table 1. The results of numerical calculations of da/dXp and their approximation in terms 
of the Peterson function are shown in Fig. 1. 

The dependence of the parameter e on the mass mi of the lightest quark can be approx- 
imated by a linear function, 

e a + hrrii (9) 

where a = 0.046 and b = 0.17 GeV^^. As the middle mass m2 is decreased, the maximum 
of the distribution da/dXp is shifted toward the largest possible relative momentum value, 
Xp = 1, while the parameter e becomes smaller. The shift of the maximum of the cross 
section da/dXp in response to the change in the mass 7713 of the heaviest quark is opposite to 
that in response to the analogous change in toi and m^' with decreasing m^, the extremal 
value of Xp decreases substantially, while the parameter e in the function given by (8) grows. 

The value e that we found here for at the mass values of mi = 0.01 GeV, 1712 = 0.5 
GeV, and ma — 1.5 GeV is very close to the values of e that were obtained in experiments 
that studied the production of Ac, Ec, and Sc baryons in electron-positron annihilation: 
e = 0.2361^1 for A^ [10], e = 0.29 it 0.06 for [11], e = 0.24 ± 0.08 for [12]. 

For baryons containing quarks of identical flavor (for example, Q^.^. ), it is difficult to take 
into account effects of interference between identical particles; therefore, the fragmentation 
mechanism cannot be so well justified theoretically in this case as for ^1^253 baryons under the 
conditions in (1). At the same time, it remains quite useful in the phenomenological aspect. 
For example, it was found in [3] that the transverse-momentum distribution of flccc baryons 
that is obtained from direct numerical calculations cannot be adequately approximated with 
the aid of the Peterson function. However, the shape of this distribution can be faithfully 
reproduced by using the so-called Lund function [13] 

D{z) ~ -(1 - zYexp{-c/z), (10) 

z 

where the parameters are set to a = 2.4 and c = 0.70. 



4. Differential cross sections and Reggeon-type fragmentation func- 
tion 

Let us now consider the representation of numerical results for the production of heavy- 
quark-containing baryons in electron-positron collisions in terms of the Reggeon-type frag- 
mentation function [14] 

D(z) ~ z^{l - zy. (11) 
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This function was employed in [15, 16] in discussing experimental data on the production of 
D and B mesons in electronpositron annihilation, the parameters j3 and 7 not being related 
there to Regge trajectories. 

An approximation of the numerical results obtained here for the production of Vlscb and 
Vlccc baryons in electron-positron annihilation with the aid of the fragmentation function in 
(11) was performed for the differential cross sections da/dxE and is displayed in Fig. 2. 

It is of interest to compare the values obtained for the parameters of the function in (11) 
by approximating the numerical results of our perturbative approach to determining cross 
sections and the values deduced from the expressions for these parameters in terms of the 
intercepts of the trajectories of appropriate hadrons. If the quark i fragments into a baryon 
consisting of the quarks i, j, and then [17] 

/3 = 1 - ftji, 7 = aqq - 2ajkq, (12) 

where the symbol q in the subscripts stands for a light quark u or d, while the quantities aij 
and ajkq are the intercepts of the meson (of quark content) and baryon (of ajkq quark 
content) trajectories. An approximate linear relation between the intercepts of trajectories 
corresponding to different quark contents was obtained within the model of quark-gluon 
strings [18]: 

2{aijk - Oiiji) = a^k - ail- (13) 

On the basis of relation (13) and the standard notation = <^dd = ctp, agg = 
Oicc — c^V" ^bb — <^T, and auud — <^udd — <^N the equalities in (12) can be recast into the form 

(3 — 1 — ax, 7 = 3(3;^ — — — 2aN, (14) 

for Qscb baryons and into the form 

P — 1 — a^, 7 = Sap — 2a^ — 2aN, (15) 

for flccc baryons. 

For some intercepts, we will take the quite reliably established values of aN — —0.4, 
ap — 0.5, and a,f> ^ [17], while, for the others, we will use the following estimates: 
a^ — —2.2 [19] and ax = —8.0 [20]. These estimates differ only shghtly from those that 
were previously obtained in [7, 21]. 

Substituting the above values into relation (14) for flscb baryons, we obtain /3 = 9.0 
7 = 4.5, but these results are in a glaring contradiction with the values of /3 = 3.3 and 
7 = 1.48 if rus = 500 MeV and with the values of /3 = 3.2 and 7 = 1.18 if = 300, 
which we obtained from a comparison of the fragmentation function in (11) with the results 
of direct numerical calculations of the cross section da/dxE- In turn, relation (15) for the 
production of VLccc baryons yields /5 = 3.2 and 7 = 6.7, which can be thought to be in very 
rough agreement with the values /3 = 2.6 and 7 = 4.6 resulting from the approximation of 
our numerical results for the cross section da/dxE- 

This suggests the simple conclusion that the perturbative and the nonperturbative (Reg- 
geon) approach do not reduce to each other and, depending on the process being considered, 
they can lead either to close or to strongly different results. 
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Fig. 1. Differential cross section da/dxp for the production of 919293 baryons in electron-positron collisions: 

(crosses) results of Monte Carlo calculations along with the errors in them and (solid curves) results of the 
calculation by formula (4) with the Peterson fragmentation bmction (8). The values of the masses of the 
quarks 91, 92, and 93 (in GeV) and of the parameter e in the Peterson function in Figs, la - 1/ are the 
following: (a) toi = 0.5 GeV, m2 = 1.5 GeV, mg = 4.8 GeV, e = 0.124; (6) mi = 0.3 GeV, TO2 = 1.5 GeV, 
ms = 4.8 GeV, e = 0.098; (c) mi = 0.075 GeV, m2 = 1.5 GeV, ms = 4.8 GeV, e = 0.066; {d) mi = 0.01 
GeV, m2 = 1.5 GeV, ms = 4.8 GeV, e = 0.048; (e) mi = 0.01 GeV, m2 = 0.5 GeV, ms = 4.8 GeV, 
£ = 0.016; if) mi = 0.01 GeV, m2 = 0.5 GeV, ms = 1.5 GeV, e = 0.24. 
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Fig. 2 Differential cross section da/dxE for the production of Q.acb baryons [(a) rrig = 0.5 GeV and (6) 
nis = 0.3 GeV] and ficcc baryons (c) in electron-positron annihilation. The crosses show the results of Monte 
Carlo calculations and their errors. The solid curves correspond to the approximation of the cross sections 
with the aid of a Reggeon-type fragmentation function (11) with the parameters (a) (3 = 3.3 and 7 = 1.48, 
(&) 13 = 3.2 and 7 = 1.18, and (c) (3 = 2.6 and 7 = 4.6. 
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